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We report a new measurement of electron antineutrino disappearance using the fully-constructed Daya Bay 
Reactor Neutrino Experiment. The final two of eight antineutrino detectors were installed in the summer 
of 2012. Including the 404 days of data collected from October 2012 to November 2013 resulted in a total 
exposure of 6.9x 10 5 GWth-ton-days, a 3.6 times increase over our previous results. Improvements in energy 
calibration limited variations between detectors to 0.2%. Removal of six 241 Am- 13 C radioactive calibration 
sources reduced the background by a factor of two for the detectors in the experimental hall furthest from the 
reactors. Direct prediction of the antineutrino signal in the far detectors based on the measurements in the near 
detectors explicitly minimized the dependence of the measurement on models of reactor antineutrino emission. 
The uncertainties in our estimates of sin 2 2#i3 and | Am 2 e | were halved as a result of these improvements. 
Analysis of the relative antineutrino rates and energy spectra between detectors gave sin 2 2#i3 = 0.084± 0.005 
and | Am 2 e | = (2.42 ± 0.11) x 10 -3 eV 2 in the three-neutrino framework. 

PACS numbers: 14.60.Pq, 29.40.Mc, 28.50.Hw, 13.15.H-g 
Keywords: neutrino oscillation, neutrino mixing, reactor, Daya Bay 


Neutrino flavor oscillation due to the mixing angle 6*13 
has been observed using reactor antineutrinos nun and 
accelerator neutrinos El 0. The Daya Bay experiment 
previously reported the discovery of a non-zero value 
of sin 2 2013 by observing the disappearance of reactor 
antineutrinos over kilometer distances EI5IEZ), and the first 
measurement of the effective mass splitting |A?n 2 e | ® via 
the distortion of the V e energy spectrum Here we 

present new results with significant improvements in energy 
calibration and background reduction. Installation of the final 
two detectors and a tripling of operation time provided a 
total exposure of 6.9xl0 5 GW t h-fon-days, 3.6 times more 
than reported in our previous publication 0- With these 
improvements the precision of sin 2 2 fl | 3 was enhanced by a 
factor of two compared to the world’s previous best estimate. 
The precision of |Ato 2 6 | was equally enhanced, and is 
now competitive with the precision of |Am§ 2 | measured via 
accelerator neutrino disappearance Gsmu. 

The Daya Bay experiment started collecting data on 24 
December 2011 with six antineutrino detectors (ADs) located 
in three underground experimental halls (EHs). Three ADs 
were positioned in two near halls at short distances from 
six nuclear reactor cores, two ADs in EH1 and one in 
EH2, and three ADs were positioned in the far hall, EH3. 
Data taking was paused on 28 July 2012 while two new 
ADs were installed, one in EH2 and the other in EH3. 
During the installation, a broad set of calibration sources 
were deployed into the two ADs of EH1 using automated 
calibration units (ED and a manual calibration system El- 


Operation of the full experiment with all eight ADs started on 
19 October 2012. This Letter presents results based on 404 
days of data acquired in the 8 -AD period combined with all 
217 days of data acquired in the 6 -AD period. A blind analysis 
strategy was implemented by concealing the baselines and 
target masses of the two new ADs, as well as the operational 
data of all reactor cores for the new data period. 

Each of the three Daya Bay experimental halls hosts 
functionally identical ADs inside a muon detector system. 
The latter consists of a two-zone pure water Cherenkov 
detector, referred to as the inner and outer water shields 
(IWS and OWS), covered on top by an array of resistive 
plate chambers (RPCs). Each AD consists of three 
nested cylindrical vessels. The inner vessel is filled with 
0.1% gadolinium-doped liquid scintillator (Gd-LS), which 
constitutes the primary antineutrino target. The vessel 
surrounding the target is filled with undoped LS, increasing 
the efficiency of detecting gamma rays produced in the target. 
The outermost vessel is filled with mineral oil. A total of 192 
20-cm photomultiplier tubes (PMTs) are radially positioned 
in the mineral-oil region of each AD. Further details on the 
experimental setup are contained in Refs. Reactor 

antineutrinos are detected via the inverse /3-decay (IBD) 
reaction, V e + p —> e + + n. The gamma rays (totalling 
~ 8 MeV) generated from the neutron capture on Gd with 
a mean capture time of ~30 ps form a delayed signal and 
enable powerful background suppression. The light from 
the e + gives an estimate of the incident V e energy, k, 
E p +E n +0.78 MeV, where E p is the prompt energy including 
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the positron kinetic and annihilation energy, and E n is the 
average neutron recoil energy (~ 10 keV). 


o Neutron from muon spallation 
• Neutron from IBD 
o Neutron from Am-C source 


a Alpha from natural radioactivity 
■ Gamma from calibration source 
♦ Gamma from natural radioactivity 
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FIG. 1. Comparison of the reconstructed energy between 
antineutrino detectors for a variety of calibration references. I?ad 
is the reconstructed energy determined using each AD, and ( E) is 
the 8-detector average. Error bars are statistical only, and systematic 
variations between detectors for all calibration references were 
< 0.2%. The ~8 MeV n-Gd capture gamma peaks from Am-C 
sources were used to define the energy scale of each detector, and 
hence show zero deviation. 

Differences in energy response between detectors directly 
impacted the estimation of Am 2 J. PMT gains were 
calibrated continuously using uncorrelated single electrons 
emitted by the photocathode. The signals of 0.3% of the 
PMTs were discarded due to abnormal hit rates or charge 
distributions. The detector energy scale was calibrated using 
Am-C neutron sources lfl8l deployed at the detector center, 
with the ~8 MeV peaks from neutrons captured on Gd aligned 
across all eight detectors. The time variation and the position 
dependence of the energy scale was corrected using the 
2.506 MeV gamma-ray peak from 60 Co calibration sources. 
The reconstructed energies of various calibration reference 
points in different ADs are compared in Fig. [T| The spatial 
distribution of each calibration reference varies, incorporating 
deviations in spatial response between detectors. Figure 1 
presents measurements of 68 Ge, 60 Co and Am-C calibration 
sources when placed at the center of each detector. Neutrons 
from IBD and muon spallation that were captured on 
gadolinium, were distributed nearly uniformly throughout the 
Gd-LS region. Those neutrons that were captured on 1 H, 
intrinsic a particles from polonium and radon decays, and 
gammas from 40 K and 208 T1 decays, were distributed inside 
and outside of the target volume. All of these events were 
selected within the Gd-LS region based on their reconstructed 
vertices. The uncorrelated relative uncertainty of the energy 



FIG. 2. Estimated energy response of the detectors to positrons, 
including both kinetic and annihilation gamma energy (red solid 
curve). The prominent nonlinearity below 4 MeV was attributed 
to scintillator light yield (from ionization quenching and Cherenkov 
light production) and the charge response of the electronics. Gamma 
rays from both deployed and intrinsic sources as well as spallation 
12 B /3 decay determined the model, and provided an envelope of 
curves consistent with the data within a 68.3% C.L. (grey band). 
An independent estimate using the beta+gamma energy spectra from 
212 Bi, 214 Bi, 208 T1, as well as the 53-MeV edge in the Michel 
electron spectrum gave a similar result (blue dashed line), albeit with 
larger systematic uncertainties. 


scale is thus determined to be 0.2%. This reduction of 
43% compared to the previous publication El was enabled 
by improvements in the correction of position and time 
dependence, and enhanced the precision of Am 2 f , by 9%. 
The reduction was confirmed by an alternative method which 
used the n-Gd capture of muon-induced spallation neutrons to 
calibrate the scale, time dependence, and spatial dependence 
of the detector energy response. 

Nonlinearity in the energy response of an AD originated 
from two dominant sources: particle-dependent nonlinear 
light yield of the scintillator and charge-dependent nonlin¬ 
earity in the PMT readout electronics. Each effect was 
at the level of 10%. We constructed a semi-empirical 
model that predicted the reconstructed energy for a particle 
assuming a specific energy deposited in the scintillator. 
The model contained four parameters: Birks’ constant, the 
relative contribution to the total light yield from Cherenkov 
radiation, and the amplitude and scale of an exponential 
correction describing the non-linear electronics response. 
This exponential form of the electronics response was 
motivated by MC and confirmed with an independent FADC 
measurement. 

The nominal parameter values were obtained from an 
unconstrained y 2 -fit to various AD calibration datasets, 
comprising twelve gamma lines from both deployed and 
naturally occurring sources as well as the continuous 13- 
decay spectrum of 12 B produced by muon spallation inside 
the Gd-LS volumes. The nominal positron response derived 
from the best fit parameters is shown in Fig. [2] The 
depicted uncertainty band represents other response functions 
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consistent with the fitted calibration data within a 68.3% 
C.L. This x 2 -based approach to obtain the energy response 
resulted in < 1% uncertainties of the absolute energy scale 
above 2 MeV. The uncertainties of the positron response were 
validated using the 53 MeV cutoff in the Michel electron 
spectrum from muon decay at rest and the continuous /3+7 
spectra from natural bismuth and thallium decays. These 
improvements added confidence in the characterization of 
the absolute energy response of the detectors, although 
they resulted in negligible changes to the measured mixing 
parameters. 

IBD candidates were selected using the same criteria 
discussed in Ref. ffl. Noise introduced by PMT light emission 
in the voltage divider, called flashing, was efficiently removed 
using the techniques of Ref. J6|. We required 0.7 MeV 
< E p < 12.0 MeV, 6.0 MeV < E d < 12.0 MeV, and 1 
/is < At < 200 ps, where E ( \ is the delayed energy and 
A/ = t,\ - t p was the time difference between the prompt and 
delayed signals. In order to suppress cosmogenic products, 
candidates were rejected if their delayed signal occurred (i) 
within a (—2 p s, 600 p s) time-window with respect to an IWS 
or OWS trigger with a PMT multiplicity > 12, (ii) within 
a (—2 ps, 1000 p s) time-window with respect to triggers in 
the same AD with reconstructed energy > 20 MeV, or (iii) 
within a (—2 /us, 1 s) time-window with respect to triggers in 
the same AD with reconstructed energy > 2.5 GeV. To select 
only definite signal pairs, we required the signal to have a 
multiplicity of 2: no other > 0.7 MeV signal occurred within 
a (i p — 200 ps, id + 200 ps) time-window. 

Estimates for the five major sources of background 
for the new data sample are improved with respect to 
Ref. gl. The background produced by the three Am- 
C neutron sources inside the automated calibration units 
contributed significantly to the total systematic uncertainty of 
the correlated backgrounds in the 6-AD period. Because of 
this, two of the three Am-C sources in each AD in EH3 were 
removed during the 2012 summer installation period. As a 
result, the average correlated Am-C background rate in the 
far hall decreased by a factor of 4 in the 8-AD period. As 
in previous publications mm, this rate was determined by 
monitoring the single neutron production rate from the Am- 
C sources. Removal of these Am-C sources had negligible 
consequences for our calibration. 

Energetic, or fast, neutrons of cosmogenic origin produced 
a correlated background for this study. Relaxing the prompt- 
energy selection to (0.7-100) MeV revealed the fast-neutron 
background spectrum above 12 MeV. Previously we deduced 
the rate and spectrum of this background using a linear 
extrapolation into the IBD prompt signal region. Here we used 
a background-enhanced dataset to improve the estimate. We 
found 6043 fast neutron candidates with prompt energy from 
0.7 to 100 MeV in the 200 ps following cosmogenic signals 
only detected by the OWS or RPC. The energy spectrum of 
these veto-tagged signals was consistent with the spectrum 
of IBD-like candidate signals above 12 MeV, and was used 
to estimate the rate and energy spectrum for the fast neutron 


background from 0.7 to 12 MeV. The systematic uncertainty 
was estimated from the difference between this new analysis 
and the extrapolation method previously employed, and was 
determined to be half of the estimate reported in Ref ©. 

The methods used in Refs. IH S) to estimate the 
backgrounds from the uncorrelated prompt-delayed pairs (i.e. 
accidentals), the correlated /3-n decays from cosmogenic 9 Li 
and 8 He, and the 13 C(a,n) 16 0 reaction, were extended to the 
current 6+8 AD data sample. The decrease in the single¬ 
neutron rate from the Am-C sources reduced the average rate 
of accidentals in the far hall by a factor of 2.7. As a result, 
the total backgrounds amount to about 3% (2%) of the IBD 
candidate sample in the far (near) hall(s). The systematic 
uncertainties in the l:! C(o,n) IG 0 cross section and in the 
transportation of the a particles were reassessed through a 
comparison of experimental results and simulation packages, 
respectively ESI. The estimation of 9 Li/ 8 He now dominated 
the background uncertainty in both the near and far halls. 
The estimated signal and background rates, as well as the 
efficiencies of the muon veto, e ;i , and multiplicity selection, 
e m , are summarized in Table [I] 

A detailed treatment of the absolute and relative efficiencies 
using the first six ADs was reported in Refs. 0121 ■ The 
uncertainties of the absolute efficiencies are correlated among 
the ADs and thus play a negligible role in the relative 
measurement of Z7 e disappearance. The performance of 
the two new ADs was found to be consistent with the 
other detectors. Estimates of two prominent uncorrelated 
uncertainties, the delayed-energy selection efficiency and the 
fraction of neutrons captured on Gd, were confirmed for 
all eight ADs using improved energy reconstruction and 
increased statistics. 

Oscillation was measured using the L/E -dependent disap¬ 
pearance of V e , as given by the survival probability 


, dn ■ i 1.267A?7i2i E 

P = 1 — cos 4 6*13 sin 2 2012 sin 2 - 


E 


, 2 . 2 1.267Am 2 e i 

- sin 2013 sin 2 - ee . 

jC/ 


(1) 


Here E is the energy in MeV of the V e , L is the distance in 
meters from its production point, 0 i 2 is the solar mixing angle, 
and A nri^i = m\ — m\ is the mass-squared difference of the 
first two neutrino mass eigenstates in eV 2 . 

Recent precise measurements of the IBD positron energy 
spectrum disagree with models of reactor V e emission BUSS¬ 
ED. The characteristics of the signals in this energy range 
are consistent with reactor antineutrino emission, and disfavor 
background or detector response as possible origins for the 
discrepancy. A separate manuscript, in preparation, will 
present the evidence in detail and provide the necessary 
data to allow detailed comparison of our measurement with 
existing and future models. Given these discrepancies 
between measurements and models, here we present a 
technique for predicting the signal in the far hall based 
on measurements obtained in the near halls, with minimal 
dependence on models of the reactor antineutrinos. In our 
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EH1 

ADI AD2 

EH2 

AD3 AD8 

AD4 

EH3 

AD5 AD6 

AD7 

IBD candidates 

304459 

309354 

287098 

190046 

40956 

41203 

40677 

27419 

DAQ live time(days) 

565.436 

565.436 

568.03 

378.407 

562.451 

562.451 

562.451 

372.685 


0.8248 

0.8218 

0.8575 

0.8577 

0.9811 

0.9811 

0.9808 

0.9811 

£m 

0.9744 

0.9748 

0.9758 

0.9756 

0.9756 

0.9754 

0.9751 

0.9758 

Accidentals(per day) 

8.92 ± 0.09 

8.94 ±0.09 

6.76 ± 0.07 

6.86 ±0.07 

1.70 ±0.02 

1.59 ±0.02 

1.57 ±0.02 

1.26 ±0.01 

Fast neutronlper AD per day) 

0.78 ±0.12 

0.54 ±0.19 


0.05 ±0.01 


9 Li/ 8 He(per AD per day) 

2.8 ± 1.5 

1.7 ±0.9 


0.27 ±0.14 


Am-C correlated 6-AD(per day) 

0.27 ±0.12 

0.25 ±0.11 

0.27 ±0.12 


0.22 ±0.10 

0.21 ±0.10 

0.21 ±0.09 


Am-C correlated 8-AD(per day) 

0.20 ± 0.09 

0.21 ±0.10 

0.18 ±0.08 

0.22 ± 0.10 

0.06 ± 0.03 

0.04 ±0.02 

0.04 ± 0.02 

0.07 ±0.03 

13 C(q, n) 16 0(per day) 

0.08 ± 0.04 

0.07 ±0.04 

0.05 ± 0.03 

0.07 ±0.04 

0.05 ± 0.03 

0.05 ±0.03 

0.05 ±0.03 

0.05 ± 0.03 

IBD rate(per day) 

657.18 ±1.94 670.14 ±1.95 

594.78 ±1.46 590.81 ± 1.66 

73.90 ±0.41 

74.49 ± 0.41 

73.58 ±0.40 75.15 ±0.49 


TABLE I. Summary of signal and backgrounds. Rates are corrected for the muon veto and multiplicity selection efficiencies e^ ■ e m . The 
measured ratio of the IBD rates in ADI and AD2 (AD3 and AD8 in the 8-AD period) was 0.981±0.004 (1.019±0.004) while the expected 
ratio was 0.982 (1.012). 


previous measurements a, model-dependence was limited 
by allowing variation of the predicted V e flux within model 
uncertainties, while the technique here provides an explicit 
demonstration of the negligible model dependence. A y 2 was 
defined as 

X 2 = - w i ■ - w t ■ N?), (2) 

where is the observed number of events after background 
subtraction in the i-th bin of reconstructed positron energy 
E rec . The superscript / (n) denotes a far (near) detector. The 
symbol V represents a covariance matrix that includes known 
systematic and statistical uncertainties. The quantity Wi is a 
weight that accounts for the differences between near and far 
measurements. For the case of a single reactor, the weight 
Wi can be simply calculated from the ratios of detector mass, 
distance to the reactor, efficiency, and antineutrino oscillation 
probability, as given by the relation: 



Here T is the number of target protons, e is the efficiency, 
and L is the distance to the reactor for a given detector. Pi 
is the oscillation probability for the v-th reconstructed energy 
bin and <f> the reactor antineutrino flux (which cancels from 
Wi). With Pi calculated in reconstructed positron energy, the 
detector response introduces small (< 0.2% above 2 MeV) 
calculable deviations from Eq. [T] 

For multiple reactor cores, the weight Wi was modified: 

»i=^ = (^)(^)E p<E ; ru ' | ri“ )r i. <4) 

The probability distribution V(Ej IUO \El ec ) accounts for the 
energy transfer from the v e to the e + and imperfections in the 
detector energy response (loss in non-active elements, non¬ 
linearity, and resolution). The extrapolation factor Tj was 


calculated as 

Tj ET es P(Ef ue , Li)<t> jk /{Liy ’ 

where P is given by Eq. [l] Lj} n ^ is the distance between 
a far (near) detector and core fc, and (j^k is the predicted 
antineutrino flux from core k for the y-th true energy bin. In 
the single-reactor core case, the antineutrino flux <j> cancels in 
the expression for r 3 and Eq. [4] reduces to Eq. [3] Although 
the cancellation is not exact for multiple cores, the impact of 
the uncertainty in reactor antineutrino flux was found to be 
< 0 . 1 %. 

The covariance matrix element V t j was the sum of a 
statistical term, calculated analytically, and a systematic term 
determined by Monte-Carlo calculation using 

% = ^ E (3 - ■ s n ($ - ■ s ?) • (6) 

Here, N is the number of simulated experiments generated 
with energy spectra S , including systematic variations of 
detector response, V e flux, and background. The choice of 
reactor antineutrino model 12214281 in calculating the covari¬ 
ance had negligible (<0.2%) impact on the determination of 
the oscillation parameters. 

Without loss of sensitivity, we summed the IBD signal 
candidates of the ADs within the same hall, accounting 
for small differences of target mass, detection efficiency, 
background and baseline. We considered the 6-AD and 8-AD 
periods separately in order to properly handle correlations in 
reactor antineutrino flux, detector exposure, and background. 
This means that i and j in the above equations ran over the 37 
reconstructed energy bins for the two near/far combinations 
and for the two periods considered (37 x 2 x 2 = 148). More 
details of this method are described in Ref. Ii29l . 

Using this method, we found sin 2 2 9 13 = 0.084 ± 0.005 
and |Am 2 e | = (2.42 ± 0.11) x 1(T 3 eV 2 , with \ 2 /NDF = 
134.6/146 (see the Supplemental Material QOl ). While we 
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use sin 2 26 12 = 0.857 ± 0.024 and A m 2 21 = (7.50 ± 
0.20) x 10 -5 eV 2 from Ref. |3T1 . our result was largely 
independent of these values. Consistent results were obtained 
when our previous methods GQ® were applied to this larger 
dataset. Under the normal (inverted) hierarchy assumption, 
|Am 2 e | yields Ato § 2 = (2.37 ± 0.11) x 10 -3 eV 2 (At 7 i § 2 = 
— (2.47 ±0.11) x 10 ~ 3 eV 2 ). This result was consistent with 
and of comparable precision to measurements obtained from 
accelerator // ;J and z/ ; , disappearance ESI HD- Using only the 
relative rates between the detectors and Ato § 2 from Ref. flOl 
we found sin 2 26*13 = 0.085 ± 0.006, with x 2 /NDF = 
1.37/3. 

The reconstructed positron energy spectrum observed in the 
far site is compared in Fig. [3] with the expectation based on 
the near-site measurements. The 68.3%, 95.5% and 99.7% 
C.L. allowed regions in the | Am 2 e |-sin 2 26*13 plane are shown 
in Fig. [4] The spectral shape from all experimental halls 
is compared in Fig. [5] to the electron antineutrino survival 
probability assuming our best estimates of the oscillation 
parameters. The total uncertainties of both sin 2 29 13 and 
|Am 2 e | are dominated by statistics. The most significant 
systematic uncertainties for sin 2 26*13 are due to the relative 
detector efficiency, reactor power, relative energy scale and 
9 Li/ 8 He background. The systematic uncertainty in | Am 2 e | is 
dominated by uncertainty in the relative energy scale. 



FIG. 3. Upper: Background-subtracted reconstructed positron 
energy spectrum observed in the far site (black points), as well as 
the expectation derived from the near sites excluding (blue line) or 
including (red line) our best estimate of oscillation. The spectra 
were efficiency-corrected and normalized to one day of livetime. 
Lower: Ratio of the spectra to the no-oscillation case. The error bars 
show the statistical uncertainty of the far site data. The shaded area 
includes the systematic and statistical uncertainties from the near site 
measurements. 

In summary, enhanced measurements of sin 2 29i3 and 



99.7% C.L. 
95.5% C.L. 



Daya Bay: 621 days j j j 

0.05 0.1 0.15 5 10 15 

Sin 2 20 13 Ax 2 


FIG.4. Regions in the Am 2 e |-sin 2 26*13 plane allowed at the 68.3%, 
95.5% and 99.7% confidence levels by the near-far comparison of 
u e rate and energy spectra. The best estimates were sin 2 2013 = 
0.084 ± 0.005 and |Am 2 e | = (2.42 ± 0.11) x 10 -3 eV 2 (black 
point). The adjoining panels show the dependence of A \ 2 on 
sin 2 2013 (top) and |Am 2 e | (right). The |Am 2 e | allowed region 
(shaded band, 68.3% C.L.) was consistent with measurements of 
| Arri 32 1 using muon disappearance by the MINOS j 101 and T2K 1111 
experiments, converted to |Am 2 e | assuming the normal (solid) and 
inverted (dashed) mass hierarchy. 



L eff /<E v > [km/MeV] 


FIG. 5. Electron antineutrino survival probability versus effective 
propagation distance L e ft divided by the average antineutrino energy 
{E u ). The data points represent the ratios of the observed 
antineutrino spectra to the expectation assuming no oscillation. The 
solid line represents the expectation using the best estimates of 
sin 2 2013 and |Am 2 e |. The error bars are statistical only. {E v ) 
was calculated for each bin using the estimated detector response, 
and L e ff was obtained by equating the actual flux to an effective 
antineutrino flux using a single baseline. 

|Am 2 e | have been obtained by studying the energy- 
dependent disappearance of the electron antineutrino inter¬ 
actions recorded in a 6.9xlO 5 GW^-ton-days exposure. 
Improvements in calibration, background estimation, as well 
as increased statistics allow this study to provide the most 
precise estimates to date of the neutrino mass and mixing 
parameters |Am 2 J and sin 2 26 * 13 . 

Daya Bay is supported in part by the Ministry of Science 
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Appendix: Why Am?.,, is used by Daya Bay 

This section describes the advantages of reporting the Daya 
Bay measurement of electron antineutrino disappearance in 
terms of an effective mass-squared difference Am 2 e , which is 
independent of the unknown ordering of neutrino masses and 
future improvements in our knowledge of the solar oscillation 
parameters. 


Introduction 

In the three-flavor framework, the survival probability of 
electron antineutrino is given by 

P(V e -A V e ) 

= 1 — cos 4 0i3 sin 2 2012 sin 2 A 2 i 
- sin 2 20i 3 (cos 2 0 i2 sin 2 A 3i + sin 2 0 i2 sin 2 A 32 ), (7) 

where A^ = Am 2 ^. The three mass-squared differences 
are subject to the constraint |Am§!| = |Am| 2 | ± |Am 21 | 
where “+”(“—”) is for the normal(inverted) mass ordering 
(or hierarchy). Therefore, determination of Am§ 2 (or 
Antjj) depends on knowledge of the mass ordering and solar 
oscillation parameters. 

The Daya Bay experiment reports a precise measurement 
of the effective mass splitting Am 2 e , which is independent of 


our knowledge of the ordering and solar parameters. In this 
approach, we approximate the survival probability using 

P(y e —> v e ) ~1 — cos 4 0i 3 sin 2 20 32 sin 2 A 23 

— sin 2 20i 3 sin 2 A ee . (8) 

Despite the advantage of using Am? e for the measurement, it 
has the disadvantage of not being a fundamental parameter. 
Therefore, we must determine a relation between A m 2 ee 
and Ato| 2 given knowledge of the mass ordering and solar 
oscillation parameters. 

In the following sections, we are going to address the 
following two questions: 

• Is Eq. [8] good enough at the current experimental 
precision? 

• How can we estimate the value of Am| once the value 
of Am 2 is obtained? 


Mathematical derivation 


Using the relation |A?n§i| = |Am| 2 | ± |Am 21 |, Eq.[7]can 
be written as, 

P(z7 e -» v e ) = 1 - 2 s 2 3 Ci 3 

+ 2s 2 3 c 2 3 y / l - 4s\ 2 c\ 2 sin 2 A 2 \ cos(2A 32 ± </>) 
— 4ci 3 Si 2 Ci 2 sin A 2 i, (9) 

where s x = sii^^, c x = cos 0 X , and <f> = 

arctan (— r ,. m2 ^ 2 U„ V The last term of the above 
\ cos 2A2i+tan^ 0i2 J 

formula is the so-called “solar term” that governs the reactor 
antineutrino oscillation at 0(100) km. For the L/E range 
covered by Daya Bay, 4 s\ 2 c\ 2 sin 2 A 2i <C 1. Thus, Eq. 
can be approximated as. 


P{v e -A v e ) 


~ 1 - 4 si.C 


13 L 13 


1 - cos(2A 32 ± 4>) 


— (solar term ,) 


= 1 — sin 2 2013 sin 2 (A 32 ± </>/2) — (solar term). (10) 


By comparing Eq.[l0|with Eq. [8] we obtain the expression 
relating Am 2 e to Am§7 (or Am^{) 

|Am 2 J = |A?n 32 | ± Am|/2 (11) 

= |Ato 31 | T (|Am 2 i| — Ato^/2 ), (12) 

where AmJ = (f> x 4P. 


Numerical evaluation 

By definition, Ais a function of L/E. Using the current 
values of Ato 2 i = 7.50 x 10 -5 eV 2 and sin 2 20 32 = 
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0.857 PTTl . Fig. |6] shows the value of Am^/2 as a function 
of energy for L = 1.6 km. We find that A.m^/2 ~ 5.17 x 
10” 5 eV 2 is essentially a constant in our L/E region, and 
numerically identical to cos^^Am^. Thus, this definition 
of A m 2 ee is similar to the definition introduced in Ref. l32l : 

A?n 2 ff | e = cos 2 0 12 \Am 2 31 \ + sin 2 0i 2 |Am 2 2 | (13) 

= |Ato 2 2 | ± cos 2 0 i 2 Ato 21 . (14) 

Figure|7]is a comparison of the approximated formula with 
Am|/2 = 5.17 x 10 -5 eV 2 , 


P ee ~ 1 — sin 2 2013 sin 2 
— (solar term), 


2 32 + 5.17 x 1CT 5 


to the three-flavor formula, Eq. [7] In this comparison, 
L = 1.6 km, sin 2 2 0 13 = 0.09, Am| 2 = 2.44 x 10" 3 eV 2 , 
and normal mass hierarchy are the inputs. The agreement 
between the two, better than 10 -4 , is excellent and exceeds 
the achievable experimental precision. 



FIG. 6 . Values of Am^/2 = |Am 2 e — Am| 2 | (black solid line) 
at L = 1.6 km as a function of the neutrino energy, with Am 2 i = 
7.50 x 10 -5 eV 2 and sin 2 2012 = 0.857 l3ll . For comparison, 
calculations based on other definitions of Am 2 e , Am 2 e = 
cos 2 0 i 2 Am 3 ! + sin 2 0 i 2 Am § 2 (red dashed line) and Am 2 e = 
4^ arcsin |^\/cos 2 $12 sin 2 A 31 + sin 2 $12 sin 2 A 3 2 j (blue dotted 
line) are also shown. 


This study demonstrates that, once we obtain the value of 
|Am 2 e | using Eq. [8] we can reliably deduce the values of 
|Ato 2 2 | and jAm^l using Eqs. 


11 


and 


12 


with 


Arn^/2 ~ cos 2 0i 2 Aro 21 . 


(16) 


Using the current values of 0i 2 and Am 21 , Am^/2 ~ 5.17 x 
10~ 5 eV 2 , and (\Am 2 21 \ - Am|/2) ~ 2.33 x 10" 5 eV 2 . 

It is important to point out that the exact solu¬ 
tion of sin 2 (Am 2 e jg) = cos 2 0 12 sin 2 (Am§ 1 + 




FIG. 7. Comparison of the survival probability at L = 1 .6 km between 
the approximated formula with Am 2 e = A 1 H 32 + 5.17 x 10 -5 eV 2 
and the exact three-flavor formula (Eq.|7}. The oscillation parameters 
used in this comparison are sin 2 2013 = 0.09 and Am 2 2 = 
2.44 x 10 ~ 3 eV 2 under the normal mass hierarchy assumption. The 
top panel shows the survival probabilities calculated with the two 
formulae, and the bottom panel shows the ratio of the two. 


sin 2 0i 2 sin 2 (Am§ 2 jg) was never used to extract the value 
of Am§ 2 or Am 2 -, from the measured A m 2 e in Daya Bay. 
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